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Abstract 

We develop a technique to formulate quantum field theory on arbitrary 
network, based on different, randomly disposed sets of scattering's. We de- 
fine R-matrix of the whole network as a product of R-matrices attached to 
each of scattering nods. Then an action for a network in terms of fermionic 
fields is formulated, which allows to calculate the transition amplitudes as 
their Green functions. On so-called bubble and triangle diagrams it is shown 
that the method produces the same results as the one which uses the gen- 
eralized star product. The approach allows to extend network models by 
including multiparticle interactions at the scattering nods. 
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1 Introduction 



Special and growing attention is brought today to quantum graphs, that 
is networks of one dimension and thin quantum wires connected at nodes. 
If such an interest can easily be explained by the current developments in 
nano-scale technology, one may also note that this domain of research is 
particularly attractive as well for condensed matter physicists as for math- 
ematical graph theory @. 

Different aspects of quantum field theory (QFT) have to be taken into 
account and developed for such a context of problem, with among others 
the behavior of fields at the junctions (vertices) of different edges. Up to 
now, it is the case of star graphs which has received special consideration 
and for which some significant results have been obtained. Such graphs are 
constituted by only one vertex connecting several edges, and can be seen 
as the fundamental parts of generic graphs. The case of a junction with 
three quantum wires has in particular been considered in the study of a 
one-dimensional electron gas [2j [3j HJ [5] with the technique of bosonisation 
naturally playing an essential role. 

Meantime, the formulation of quantum scattering theory on a graph, 
based on the connection of self adjoint extensions of the Schrodinger op- 
erator and Hermitian symplectic forms, provided an explicit expression of 
the unitary scattering S-matrix formed by the reflection and transmission 
amplitudes and encoding the allowed boundary conditions at the vertices 
[6]. Moreover, in the works [BJ [TJ [8] different, but equivalent rules for ob- 
taining the transmission and reflection amplitudes were given. One rule is 
based on the so called generalized star product. All rules are based on linear 
algebra and are relatively simple to apply. In particular they may easily be 
implemented for the purpose of numerical calculations. 

Then, combining these last results with an algebraic field theoretical 
treatment for quantum integrable systems in (1+1) dimensions with a re- 
flecting and transmitting impurity [9] allowed the authors of \10\ [TT1 [T2] to 
develop a general framework for constructing, via the so-called reflection- 
transmission or "R-T" algebra, vertex operators from bosonic fields prop- 
agating freely in the bulk and interacting at the junction. After a de- 
tailed study of scale- invariant interactions, and an explicit determination 
of the critical points for a star graph with any number n of edges, the 
four fermion bulk interaction could then be considered and the Tomonaga- 
Luttinger/Thirring model solved at such critical points. This result is finally 
used for determining the charge and spin transport as well as to establish 
a relationship between them, recovering and extending the previous results 
first obtained in [2] , and later confirmed in [3 IH E] . 

Such interesting results urge on the study of more complex graphs rep- 

4 see for example [T] for a recent and detailed overview on this subject 
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resenting more realistic situations, the rather obvious next step being the 
case of a network with more than one vertex, linked together by internal 
wires, and particularly configurations where closed loops are present. The 
simplest examples in this last case are the "bubble diagram" (see below 
Fig. 7) with two vertices related by two internal wires, and the "triangle dia- 
gram" (see Fig. 8) where three star graphs combine in such a way that each 
couple of vertices is connected by an internal line. In order to extend the 
previous formalism developed in [10\ \TT\ [12], one can wonder whether the 
closed diagram made of several junction points, each one connected with 
an internal as well as one or more external wires, can be represented by an 
unique "mega" vertex with the same external lines. This is one of the re- 
sults of this paper, where an explicit construction of the scattering S-matrix 
relative to the "mega" vertex is given in terms of the scattering S-matrix 
elements of the different vertices involved. At this point let us mention that 
direct computations of the final S-matrices have recently been performed in 
the case of a chain of impurities Ref . [13] as well as for closed loops Ref. [H] . 
A direct comparison of some of the results obtained by this last approach 
with ours can easily be done. But we must stress that the technique we are 
using is completely different from the one of Ref.[13t HH] and will allow us, 
as explained below, to propose an action formulation for models defined on 
a lattice associated to the considered network. 

It is the extension of transfer matrix approach which is hereafter used 
to handle this problem. More precisely, to each of vertexes of the network 
we allocate a R-matrix, which is defined by the corresponding scattering 
S-matrix. Then we construct a product of all R-matrices of the network 
in an appropriate way defining mega-R-matrix. The resultant R-matrix 
defines an S-matrix of the whole network. This is one of main results of the 
paper. Then, using an approach first initiated in Ref.|15j. we rewrite the 
R-matrix constructed for each star graph in terms of creation /annihilation 
operators ( representing the in/out-coming particles) and with coefficients 
directly related with S-matrix entries. The replacement of such operators by 
Grassmann-valued fields in 2-dim. space allows one to perform directly the 
product of such R-matrices, each associated to a different vertex connected 
to the next one by (at least) one internal line, and finally to reconstruct 
from the so-obtained transfer matrix the one mimicking the scattering of 
the "mega" vertex containing the different connected junctions. 

This technique, namely the transformation of the R-matrix into its coher- 
ent states basis in fermionic Fock-space allows to give an action formulation 
of the corresponding model on a lattice formed by the network. The explicit 
check of equivalence of the action formulation of network models on a basis 
of R-matrix with the star-product rules of calculations, earlier formulated 
in the paper [TJ, for so-called bubble and triangle diagrams is the another 
important result of this article. Note that the general proof of equivalence 
is given in [16] . 
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The action formulation is a convenient way, if not the only, to investigate 
networks with large number N of vertices, or near the critical point, when one 
can develop an equivalent quantum field theory approach to the problem. 

The importance of the action formulation of network models is hard to 
overestimate. It can be applied, for example, to the Chalker-Coddington 
(CC) phenomenological model [T7j devoted to plateau-plateau transitions in 
quantum Hall effect. This approach to the CC model was developed in |18j . 

Let us at this point note that in Ref.[15] such a fermionization tech- 
nique was also used for the formulation of lattice models in (2+1) dim. in 
connection with the 3D-Ising model. Applications for the construction of 
integrable models with staggered parameters have also been found [19]; see 
also |20j for more developments. Obviously, the direct connection between 
such lattice models and quantum graphs has to be exploited more deeply, 
and we will return to this point in our conclusion. 

The paper is organized as follows. We start in Section 2 by considering 
the case of several impurities on a line. Then one easily notices that, by 
simply re-organizing the entries of the scattering S-matrix relative to each 
vertex, one gets transfer (T-) matrices such that their product, following the 
order of the different vertices along the line, provides one with a transfer 
matrix relating the external legs of the diagram: it is then a simple exercise 
to deduce, by an operation inverse to the previous one, a matrix which can 
be seen as the S-matrix of the "mega" vertex of the diagram. Actuallythe 
T-matrix associated to each vertex can be seen as the one particle sector of 
the extended XX-model R-matrix. 

In the Section 3 we give the formulation of three-channel -matrix 
in a standard matrix and operator forms. Then we pass to coherent state 
basis and present a simple expression of the R<® -matrix via three channel S- 
matrix. Then we formulate a problem of complicated networks and express 
total i?-matrix as a product of its ingredient and R^ matrices. This 
allows us also to give an action formulation for any complicated networks in 
via fermionic fields. 

In Section 4 we apply our method to calculation of the so-called bubble 
diagram, while in Section 5 we present the calculations for the triangle vertex 
diagram and show, that corresponding mega-S-matrices coincides with the 
results obtained via star-product procedure given in [61 [7] . 

2 Two particle R-matrix 

2.1 Connection of Transfer matrix with R-matrix 

The act of the scattering of the particle on a potential center /impurity can 
be described by a 2x2 S-matrix of the form 
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where i = 1,2 denote the number of channels: Figure 1 demonstrates this 
process. 



in out 



out m 

A 

Figure 1: Two channel S-matrix. 

The S-matrix consists of transmission T\ 2 (k), T 2 \{k) and reflection Rn(k), R 2 2{k) 
amplitudes ( details can be found in [U [T3J [21] . 

/ R u (k) T 12 (k) \ 
b ~ { T 2l {k) R 22 (k) ) ■ {2 - 2) 

This two particle process can be described also by the transfer matrix 
which maps the states | 1 >j n , | 1 > ou t on the left hand side of the scattering 
region A to the | 2 >„, | 2 > ou t on the right hand side. 




T 12 (k) 
fln(-fc) 
"Tia(-fc) Ti 2 (-fe) 




(2.3) 



In ()2.3p the symmetry properties (unitarity, analyticity, consistency) of 
the S-matrix (I2.2p have been taken into account. [6| I21|, fl3] . 

Since the transfer matrix maps the states on left hand side to the states 
on the right hand side of the scattering region the process of scattering on 
the multiple points will simply be described by the multiplication of the 
transfer matrices (see Fig. 2). Continuously mapping the states from the left 



Tfflinr 



Tj T 2 T 3 T 4 T n 

Figure 2: Chain of two channel T-matrices representing a chain of equidis- 
tant impurities. 

to the right on the Fig. 2, we will have a product of transfer matrices 

r (n) =T n T n _ 1 ...T 1 , (2.4) 

the matrix elements defining final transmission (k) and reflection (k) 
amplitudes. This procedure of calculation of amplitudes is much simpler 
than calculations in S-matrix formalism, where we do not have local matrix 
multiplication but should deal with counting of non-local effects. 
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It appears that this transfer matrix can be considered as an one particle 
sector of the R-matrix of an XX-model. Namely, the standard form of the 
extended XX-model's R-matrix can be represented as: 



Rxx 



where we denoted 



R oi 

E>00 



T22 














1 










rn 


rnr 2 2 - r 12 r 2 i 

















\ 



(2.5) 



-Til, 



plO 



T\\T22 ~ T12T21, 
Tl2, 



r\{ -- 

nOl 
^10 



-T22 
-T21- 



(2.6) 



Then, by use of unitary properties of the S-matrix (|2.2I) one easily can see, 
that the middle block of elements R\q, Rq\, Roi> R% wm coincide with the 
transfer matrix (|2.3jl after dividing them on —Rqq 

T 2 i(k) = m, T 12 {k)=r 22, 

Rn(k) = r 21 , R22(k)=r 12 . (2.7) 
Hence, the r-matrix elements are connected with the S-matrix ones as: 

r = C-S (2.8) 



with C 




the numeration of spaces of the initial states differing 



in the R and T matrix formulations (compare Fig.l with Fig. 3 below). 

Via Jordan- Wigner transformation, the R-matrix (|2.5p can be repre- 
sented as an operator consisting of fermionic creation/annihilation operators 
which act on a corresponding Fock space made of two quantum spaces (see 
Fig. 3 below). 



R<$ : V x <g> V 2 -> V x <g> V 2 . 



(2.9) 



Now, let us introduce scalar fermions cf, Ci in each linear space Vi, i = 
1, 2 forming a set of quantum states \i >i n / ou t- Defining the R-operator as 



R 



(2) 
12 



=: e 



(2.10) 



where we have normal ordering convention of operators in space V\ and 
anti- normal ordering convention space V2 , one will obtain [22} [23] 



R 



(2) 
12 



where n 



R\\nin 2 + -Roonin 2 + -^io n i n 2 
+ R%\ n x n 2 + R% 4 c 2 + R\\ 4 ci , 

1, 2 is the particle number operator and hi 



(2.11) 



m. 



5 



2.2 Chain of Equidistant Impurities 

The transfer matrix approach allows to analyze very easily a chain of scat- 
tering centers/impurities and to calculate the transmission/reflection am- 
plitudes. In this section we will explicit such a calculation by considering a 
chain of equidistant, identical impurities. 

For such a purpose, we start by defining the local transfer matrix Tj at 
the chain point i, with i = 1, ...,n by translating the transfer matrix To 
defined by the formula (|2.3p at the point zero via translation operator 

/ ikxi n \ 

"* = ( <> e-°-.j < 2 - 12 > 

as Tj = W~ TqW^. Then the final transfer matrix appears as the product 
of the different local transfer matrices Tf 

T (n) = T n T n _x ...T\ = W-^TqWTW-'Wu (2.13) 

where W = Wi+iW^ 1 is constant due to equidistant position of impurities 
d = Xi + i — Xi = const. 

We need to calculate T". In order to reach this goal one can diagonalize 
the matrix TqW with the help of an unitary matrix U: 

T {n) = W-'UiToWr^^U+W-'W, (2.14) 

The eigenvalues and eigenvectors of TqW can be obtained via: 



e ikd R 11 (k)e~ ikd 
T^k) T 12 (fc) \ ( U± \ _ ^ ( U± 

Ru(--k)e M °~ ikd 




A ± | ™ , (2.15) 



T 12 (-k) T 12 (-fc) 

where eigenvalues X± easily can be found with the characteristic equation: 



(— 



ikd Q—ikd 



x - x [nm + w) +1=0 - (2I6) 

Defining 



a(k) = ArcCosh 



(2.17) 



2 \T 12 (k) T 12 (-k) / 
one can represent eigenvalues in a form 

A ± = e ±a ( fc ). (2.18) 
After some more algebra, one gets for T^: 
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rp(n) 



sinh(na) 1 ik(n-l)d 
sinh(a) T^^F) 
sinh((n-l)g) c ik(n-2)d 
sinh(a) 

sinh(na) Rn(-k) ik(2xi +(n-l)d) 
' sinh(a) T 12 (-k) e 



V 



sinh(na) giiffc) -jfc(2g] +(n-l)d) \ 
sinh(a) T 12 (fc) C > 



sinh(raa) 1 -ifc(n-l)rf 
sinh(a) IW-fc) 

sinh((n-l)a) -ik(n-2)d 
sinh(a) 



(2.19) 



/ 



Note that, in the case of only one impurity, i.e. n=l and d=0, one 
directly recovers the T-matrix defined in Eq. (2.3) - as it should be. 



2.3 Fermionic Coherent States 

In this subsection we will define fermionic coherent states which are very 
instrumental in constructing of action/Lagrangian formulation of network 
models |15j . They form another basis of the linear space of the quantum 
states, giving the possibility to express the matrix elements of the R-matrix 
via Grassmann variables. 

By definition, fermionic coherent/anticoherent states are eigenstates of 
fermion creation/annihilation operators c, c + . Since the operators c, c + have 
anti-commuting properties, their eigenvalues cannot be ordinary numbers. 
They will be Grassmann variables satisfying the following properties: 

ifjlj) = Iplfj = 0, IpTp = —tpll'} 

J dtp = J dtp = 0, J d^jip = J dijtp = 1. (2.20) 

With the help of these variables, let us define a coherent state: 

|V> = e-^ c+ |0), (V5|=(0| e -^, (2.21) 

as being an eigenstate of the annihilation operator c with the eigenvalue ip: 

c | if)) = ip | if,) | c + = | ip. (2.22) 

Correspondingly, the anti-coherent state 

|$ = (c+-$|0), (V> |= (0 I (c + V), (2-23) 

is an eigenstate of the creation operator c + with the eigenvalue ip: 

c + \ ip) = ip \ ip) (V> | c = -(i/> | ip. (2.24) 

By use of the relations (|2.20[) . it is straightforward to check that we get the 
following normalization: 

$ | V) = (V 1 I V 5 ) = (2.25) 
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and the completeness relations: 

j ##|^)<$|e^ = l, Jd^dtj} | ^)(V> | e** = 1 (2.26) 

for the coherent and anti-coherent states respectively. 

The coherent states are forming a new basis in the two dimensional 
Fock space of fermions, where it becomes very easy to calculate the matrix 
elements of any functional in the operators c, c + written in a normal/anti- 
normal ordered form. In particular, it is straightforward to find out that, 
for any normal ordered operator : F(c + , c) : 

$ |: F{c + ,c) :| V) = F($,^)e^ (2.27) 

while for an anti-normal ordered operator: 

(ip \: F{c + ,c) :| -0) = F§, -tp)e^ . (2.28) 

Note also that the kernel F(tp, ip) of two operators : F{c + , c) :=: Fi(c + , c) :: 
F2(c + , c) : reads: 

F(^) = | dxdxe-™F 1 $, X )F 2 (x,il>)- (2-29) 

We will use this relation extensively below. 

The R-operator (|2,lip has a very simple form in the coherent state ba- 
sis. Introducing coherent/anti-coherent states tp[, tpi, ip2, m the spaces 
V1/V2 correspondingly defined by formulas (|2.2ip / (|2.23j) . one obtains for the 
matrix elements (ip^j^'i I R I V'liV^of the R-operator: 

(tp2,i>i\R \ ^1, ^2) = e ~ ril ^ 1 ~ r22 ^ 2 ^~ ri2 ^ ! 2~ r21 ^ 1 

= e -i>i(.cs)iji)j _ (2.30) 

In order to make a correspondence with the notations on Fig.l, where 
the left hand states are marked with index 1, while right hand states are 
marked by 2, we redefine ip\ = rji, ip\ = 772, tp 2 = ^'21 ^2 = an d obtain for 
the i?-matrix the expression 

( V '2,m I R I VuV'i) = e~^ s ^. (2.31) 

Hereafter, we shall use i?-matrices in this last formulation like (I2.31|) . 

Graphically, the R-operator in the coherent state basis can be repre- 
sented as in Fig. 3, where it becomes clear that the S-matrix parameters are 
merely hopping parameters. 

Note that, in getting formula (|2.3l)p we have used formulas (j2.27H2.25]) 
and (12TTU1) . 
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2 — 
V 2 



V 2 



Figure 3: Convenient representation of two-channel R-matrix as an quad- 
rangle where the sides are hopping links of fermions. 



3 Three particle R-matrix 

3.1 Three particle S-matrix 

In this Section we will define the three particle R-matrix connected to the 
corresponding three particle S-matrix. The main problem that we would 
like to solve is the following. Usually, in order to describe the scattering 
of three particles we define a three particle S-matrix. But then it is very 
hard to make calculations in a situation when of a network with multiple 
scattering points. The problem stands in non-local rules of the S-matrix 
multiplications. The situation is totally different if we have to deal with 
multiplication of R-matrices. Then we have local multiplication rules of the 
tensorial type allowing to handle complicated problems with networks. 

We will start from the three particle of S-matrix and define a correspond- 
ing R-matrix in the next section. 

So, let us consider three channels where the incoming and outgoing states 
are marked as | i >i n and | i > ou t, (i = 1)2,3) respectively. In a way 
similar to the two particle case (12. If) , the three particle S-matrix is an matrix 
mapping incoming states to outgoing ones: 

I i ^ > out = Sij | j >m (3.32) 

It is convenient to represent this scattering process by the vertex shown 
in Fig. 3. 



3.2 Definition of Three Particle R-matrix 

Up to this point the S-matrix was describing two channel scattering processes 
where 0,1,2 particles could participate. Now we introduce a three channel 
R-matrix containing information about 0,1,2,3 particle processes involved in 
the scattering. 
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out 




Figure 4: Three channel S-matrix/R-matrix. 

For further convenience, we change the traditional notations of the states 
in the channel presented in the Fig. 3 and make the re-numeration in accor- 
dance with Fig. 4. Namely, let us change |2') — ► and — ► |2') and 
introduce the matrix: 

r = CS, (3.33) 
where S is the three channel S-matrix (|3.32p and the matrix C 



C = 1 . (3.34) 




ensures the re-numeration of the states according to Fig. 4 and Fig. 5. 

The three particle R-matrix is a matrix/operator acting on the direct 
product of the three two dimensional linear spaces Vj, j = 1, 2, 3 ( Fig. 5) 

#123 : Vi ® V 2 ® V 3 -► V x ® V 2 ® V 3 . (3.35) 

The dimension of V spaces is defined by the number of incoming and out- 
going states. In the simplest case it is two and we will consider this case in 
this article. Let us also define the basis states for the each space Vi as |0)j 
and 

Our aim is the definition of the R-matrix corresponding to the given 
S-matrix. Any R-matrix can be represented in three different bases of the 
linear spaces Vi, i = 1,2,3. The first, and the most popular one is an 
ordinary matrix form. In case of three particle, the R-matrix entries read 
i?^^ 7 where each index takes the values and 1. As we will see below 
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this notation is justified because it is connected with the particle (fermion) 
number in operator representation of the R-matrix. We will define and 
consider R-matrices corresponding to the processes with particle number 
conservation. The most general form of this type of R-matrix was defined 
in [20] and have a following form: 



' -^000 



R 



\ 



pOOl 
-"■001 



pOOl 
-"010 



pOOl 
-"100 



pOlO pOlO pOlO 
-"■001 -"oio -"100 



plOO 
-"■001 



plOO 
-"010 



plOO 
-"100 



pOll pOll pOll 

-"■on - fl ioi -"no 



plOl 
pllO 

-"on 



plOl 
-"101 
pllO 
-"101 



plOl 

-"no 

pllO 
-"101 



V 



pill 



(3.36) 



It is presented in a block-diagonal form, where each block corresponds re- 
spectively to the 0—, 1—, 2—, 3— particle scattering: we interpret |0), |1) 
states as vacuum state and one particle state. 

Each two dimensional space Vi,i = 1,2,3 can be realized as a Fock space 
of scalar fermions. This will allow to formulate the operator R-matrix in a 
fermionic space equivalent to the Jordan- Wigner transformation. 

Let us introduce the creation/annihilation operators Cj, cf of scalar fermions 



in each linear space Vi, 



1, 2, 3 such as: 



Ci \0)i = 0, c+|0); = l^i- 



and consider the operator [20 



R 



123 



,(-W(*y+ry)«y 



(3.37) 



where the matrix elements rij are forming the set of parameters of the model 
and are connected with S-matrix entries via formula (|3.33|) while the symbol 
:: means normal ordering of the operators in the space i = 1,3 and anti- 
normal (hole) ordering for the fermions in i = 2. 

Expanding the exponent in the formula (j3.37H and after some algebra 
one can obtain the following expression for the operator R-matrix: 



^123 = #000 n l n 2™3 + -R 01 n l n 2™3 + # 10 n l n 2™3 + -Rl n l n 2n 3 

+ Rol\nin 2 n 3 + R^n^n^ + R^n^n^ + Rl\{nin 2 n 3 

+ (R^c+cz + ^ lC +c 3 )ni + (Rl$c+c 2 + RZc+c 3 ) ni 

+ + R%4c 3 )n 2 + «?4ci + R%4c 3 )n 2 

+ («o?o4 ci + R^4c 2 )h 3 + (R$l4ci + ^?oic+c 2 )n 3 (3.38) 
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with 

Rill = r n r 33 - r 13 r 31 , = r u r 2 2 - r 12 r 21 , Rool = r 33 r 22 - r 23 r 32 , 

RZ = 1, R\ \ = -Det[(-iy rij ], R\\ = -r u , JFaggg = -r 22 , R° \\ = -r 33 , 



R iw = r n r 23 - ri 3 r 2 i, R^ = r 12 r 3 i - r u r 32 , R( 01 = r V3 r 32 - r 33 r l2 , 
R oii = r 33^2i - r 23 r 31 , RiqI = r 22 r 13 - r u r 23 , = r 22 r 31 - r 2i r 32 , 

^110 - r 13, ^011— r 31, RoiO - _r 23, 

Roll = r 32 , RTo = rw, RZ = -r2i- (3.39) 

In the spaces Vi, i = 1,2,3 the creation-annihilation operators can be 
identified with |a)i(a'| as: 

* = \0)i(l\, 4 = \l)i{0\, 4a = \l)i{l\, acf = |0)i(0|. (3.40) 

Note that with this definition, each Vi space acquires a grading: \a)i\(3)j = 
(— l) PaP f\l3) j\a)i, i 7^ j, where p a = a is the parity of the state \a). 




time 



Figure 5: The R-matrix, attached to the S-matrix and defined by this pic- 
ture, has the entries R^pL 1 

In this ket-bra language, the i?i23-operator (|3.37|) can be represented via 
its matrix elements R°^^ = (jPct \ R \ a'P'^y 1 ) in the following way: 

^123 = </ 7 Y (l«)il/?)2l7)3)(3(7 / |2(/3 / |i(«'|) (3.41) 

= ( - 1 )+ ^' ( ^ } ( I «> ! I ) ( I /?> 2 I ) ( I T> 3 <T' I ) - 

The operator -R123 has even parity, and for its non-zero matrix elements the 
relation (-i^+^+^+JV+P^+Py = 1 takes place. 
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The R-operator appears to have a very simple form in a coherent space 
representation defined by formulas (|2. 21112. 23]) . By using a coherent state 
basis (|2,2ip for the spaces 14)^3, and an anti-coherent basis (I2.23P for V2 as 
well as the relations (|2.27l I2.28[) , one can obtain a very transparent connec- 
tion between the R- and the S-matrix. Namely: 

R = e -An^ = z-viSijVj j (3.42) 

where we have summation over i,j = 1,2,3 and r/j = ifti, fji = Cijtpj. 
According to the definitions in (|3.33[) . the first equation in (|3.42p corresponds 
to the Fig.5, the second one to the Fig.4, i.e and S%j are the hopping 
parameters defined on the network of Fig.5 and Fig.4 respectively. 

One needs also to define the -R123 operator with inverse orderings: hole- 
ordering and anti-coherent basis for the spaces V\ , V3 , and normal ordering 
and coherent state basis for V2. This will bring to the respective changes of 
the expressions of the matrix elements of R^ written in terms of nj matrix 
elements. 

The connection of three particle R^ -matrix (|3.33p with two particle one 
defined by (|2.6|) is very simple. Namely: 

K af3 ~ K a ( 6A6 ) 

With this three-channel R-matrix R^ and the two channel R-matrix 
i?( 2 ) defined in Section 2 at hand, we are in position to consider complicated 
networks with any number of R^ and R^ vertices. The generalization to 
multivertex (n > 3) processes is also straightforward: 

Ru...n =■ e ^ { - 1)3 ' c t {Sij+rij)c i := (3.44) 

= R a\a 2 '.'. a a n (|"l)l|"2)2 • ' ' |«n)n) (n(a' n \ • • • 2 (a' 2 1 1 |) 

Here also for further convenience one can define two ways of the orderings: 
normal ordering and hole ordering for the operator expressions respectively 
in the spaces V2k+i and V2k and the inverse situation. For the first case one 
must take f = j in (|3.44|) . for the second case f = j + 1. 

3.3 Networks 

In this subsection we will describe a procedure for handling a process with 
multiple external channels and many R^ and R^ vertices. Our aim is 
to express the i?-matrix of the whole process via its ingredients. Let us 
consider now an example of complicated network as presented on Fig. [6l 
We have attached the indexes a, (3, 7, 5 u, p, (ft, lu to the 8 incom- 
ing external states and ct\, {3±, 71, S\ v\ % p\ <j>\, uj\ to the 8 outgoing 
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external states of the network, while the internal indexes are denoted by 
X, Xi, h, fa, fa, S 5 , v 2 , u 3 , 1/4, 1/5, 0, 

61, a, <7i, pi, P3, P4, (f>2, <f>3, 04, t, T\. There are eight three particle 
and two particle R^ matrices here disposed at the points with coordinates 
£i, i = 1 • • • 10 in the space. 




Figure 6: Example of network. 

But the definition of the two- and three channel i?-matrices given above 
doesn't depend on the coordinates of the space; they are defined at the zero 
point. Therefore, in order to form a network with different positions of the 
scattering nodes, we should translate each of the scattering matrices from 
the zero point to their position at £. This can be done by use of the set of 
translational matrices W^ki), one for each of channels (i = 1,2,3), defined 
as 

with ki = krii, where n*j is unit vector exiting from the vertex £ along a 
channel i = 1,2,3 in direction to the neighbor vertex (defined by arrows 
on the figures). The matrices W^ki) and W^ki)~ l = W^(-ki) are acting 
on the incoming and outgoing states of the channels respectively. On the 
matrix elements (I3.39D of the i?-matrix they act as: 

• (^))i(^- 1 )(fe))£(W | <fc3)))i (3.46) 
In the fermionic operator formulation (|3.37p . this transformation reads as 
#123$ = (= e-V-'-^rt* :) (: e^ 1 ^ :) (: e -(l-e- e1 *)4<* .) Ru3 
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x(: e -(l-^)c+ca :) ( . e (l-e-^2) c + C2 ;) ( . e _(i_ e «V3 )c + C3 ;)) (3 _ 47) 



which is equivalent to the corresponding change of the parameters 
(ICTD . Namely, 



nj m 



r«(D = e-^-^nj, (3.48) 

where the are defined in f|3.33p . (|3.42l) . 

The equation (|3.48p gives also the transformation formula for the two- 
channel matrix elements of two-channel R^^p . 

With the notations of R^ -matrix elements (I3.39p . the resulting R^- 
matrix elements, describing the process of Fig. 6, can be written as 

(8) ai a Si "Pi P = J_ ^(iW^+Px+Pe 

7 71 S ■uji p p% (p ipl jy Z_^\ x / 

„(3)<Wi ( 3 )<W B (3)XlS 3 V3 p(3)Xai« E>(2)<^4 f . s 

(3)^1^1 „(3)¥>3P30 „( 3 )W2P2 „(2)<P4P4 ^yipri 

11 5 i U i 9 JrL ipifii9 1 Jrl VDx^pzP'i <P2P2 V4P4T 

where we have summation over all internal repeating indexes. 

The normalization factor N in the expression (|3.49|) is essential to factor 
out vacuum(closed loop) processes in the network, which are going on with- 
out inclusion of external states. This factor will be concretely specified in 
examples presented in the next two sections. The signs (— l) p are reflecting 
the graded character of the i?-operators (|3.44j) : for recovering the formula 
(I3.49p . one has to make the product of the operators (|3.44p corresponding 
to each vertex and to sum up over all the internal states. 

One must also takes into account the following multiplication rules: if 
the product of two operators relative to a given internal channel n contains 
a product such as n {a\ x \a) n , it must be evaluated using n (o:\oi) n — 5 Q q, and 
if a product like \(3) n x n (j3\ occurs, then it can be derived using: \f3) n (f3\ = 
(_ l)P/3(5g^. This second relation can be obtained formally from |/^)^(/3| — 
(— l) P/3 n (/3 |/9) n , or exactly by introducing a supertrace over the states in the 
n-th channel. 

It is clear that this procedure to express the total i?-matrix via its con- 
stituents can be generalized to any complicated network. Below we are going 
to apply this technique to the bubble and triangle diagrams. 



3.4 Action Formulation of the Networks 

In this section we formulate the network problem in terms of the integral ker- 
nels of the -R-operators ( formulas 1 2 . 3 1 3 .4*21) in the space of the Grassmanian 
functions and express an action principle for it. 

As in subsection 3.2, we attach to each local vertex with m-channels a 
fermionic vertex-function: 

Rmd^iM) = e-Xr^i, (3.50) 
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which is the kernel of the corresponding .R-operator (|3.44p . the being the 
entries of the rmatrix describing the scattering process over that vertex. In 
order to find the r—matrix(5'— matrix) corresponding to the entire network, 
one must take the product of the -R-operators (|2.3QI3.42p corresponding to 
the vertexes and integrate in the functional space over the internal fermionic 
fields. 

Suppose that we have a network with n vertices, each of which with 
Si, i = 1, ... n internal links (which connect two vertices of the network) 
and qi external links (which are ending or starting in one vertex). Any link 
in the diagrams has its double with an opposite direction for the arrow). 
Let us denote the Grassmann- variables connected with the links of the k-th. 
vertex as ipf, ifjj: by a convention the variables will be attached to the 
links with outgoing arrows, while the variables ipj attached to the links with 
incoming arrows. Every internal link connects two vertices possessing two 
pairs of the Grassmann- variables: in the diagram below (Fig. 7) they are 
rff, ip!? and V$', ipf, . 




Figure 7: Attaching of the Grassmann-variables to the network. 

Note that we are considering general networks, i.e. without any restriction 
on their structures, and, therefore, any two vertices can be connected by an 
arbitrary number of links. Below we denote the set of Fermi-fields on the 
internal links by Ui n , and the set of Fermi-fields on the external links by U ex . 

Taking into account the integration rule for the product of integral ker- 
nels ()2.29p . the vertex- function associated to the network will be: 

, Y>f,'} € U ex ) = e~ £m/ S.^fe $i s »> ^ = 

i „ n n 

= n 11 U d ^ d ^I[ R ^t4}) (3-51) 

{i,i'}£U in k,k' k 
N {i,i'}eU in k,k' 

Here Ui n denotes also the set of internal links between the vertices with 
dim\Ui n ] = Y^j s j = 2s, where s is the total number of the internal links in 



16 



the network. Similarly U ex is the set of external links and has dimension 
q = YskQki that is the number of the external links of the network. The 

indices k, i and k' , i' in the measure of integration e~^ i correspond to the 
end points of the same link in the network (hence the mark ' appearing on 
the symbols in the products of (I3.51|) and in the sums of (I3.52|) below). The 
normalization factor N is present in order for the mega-network i?-matrix 
to appear in an exponential form (|3.5U|) with an identity pre-factor. It is 
defined by the contribution of the closed loops to the functional integral 
(13511) . 

The argument of the exponential function in the integrant in (I3.5ip : 

' n n 

a{{$, fitfy g u m ffi u ex ) = - E E ®$ - E E rWflij 

i,i'&U in k,k'=l i,jeUex(BUin k=l 

(3.52) 

can be regarded as an action, and the normalization factor N as a partition 
function for the corresponding network model with fixed external (bound- 
ary) fields. This interpretation allows also to calculate the 5-matrix elements 
of the total network as a Green function of the external Fermi-fields: 

{Q)% = (S- 1 )™: (3.53) 

Now we will show that presented rules of calculation of the S-matrix 
of the mega network reproduces the star-product rules formulated in the 
articles [7]. A general proof of equivalence of two ways of calculations is 
available in a paper [16]. Here we will demonstrate their equivalence for 
so called bubble and triangle diagrams. Moreover, we show that the same 
result can be obtained via direct R-matrix products by following the rules 
presented in the example (|3.49[> . 

Let us separate the summands in the exponent which contain variables 
with internal indexes and represent the integral (|3.51|) in this way: 

i - e,, 3£Wci ei / n n ( 3 - 54 ) 

{i,i'}&U m k,k> 

(k) 

where A r is a 2s x 2s matrix, with the elements A rik j k , = Si k j k , + 5kk'r\j , 
where in the second summand only the elements relative to the internal links 
are taken into account. This matrix, which contains only unit elements on 
the diagonal, can be presented in a block form after some rearrangement 
of the rows and columns , with n diagonal blocks containing the parts of 
the matrices, referred to the internal links of the network, and 2s unity 
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elements disposed out of the diagonal blocks in such a way that each column 
and each row contain only one unity element. When "tadpole" -vertices are 
present in the networks, elements of the corresponding matrices also 
appear in the diagonal together with the unity elements. 

After integration over internal tp' 's in the Gaussian integral (I3.54p . we 
arrive at: 

N 

from which it follows that the normalization factor N in (I3.51|) must be 
equal to det A r and the mega-S'-matrix is: 

Sf = r^S kk '-A r f, (3.56) 

where 



and k, k' are relative only to vertices which possess external legs. 

Expression (I3.55I) easily follows after using Hubbard-Stratanovich trans- 
formation in (|3.54p and: 

2s / q n \ ( q n 

i'j'eUin yi&u^ k J \j£U ex k' 

q n 

= E (3-58) 

i,j&Uex k,k' 

For the particular case n = 2 let us write down a detailed expression of 
the network vertex-function (or the kernel of the i?-operator corresponding 
to that network). Now the network contains two vertices having respectively 
s±, q\ and S2, q2 internal and external links. After dividing each of the 
matrices into four blocks, following the definitions in the work [7], the kernel 
function of k-th vertex can be represented as: 

g ij r i r ] — g ij 1 3 

Here VE^ = {ifii, • • • , V's,.} 1 " is the set of variables connected with the internal 
links of the fc-th vertex, and ^ = {ipg k +ii " " i V's fc +g fc } r the set of variables 
connected with the external links. Correspondingly: ^f' = {y>i, • • • , rpg }, 

/ .(*) .(*) \ 

*2 = ■ ■ ■ » iPs h + qk }- I* follows that: r« = & g . 
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Now the expression (I3.51|) takes a simple form (it is obvious that si = 
s 2 = s): 

f d^\d^\d^ld^\ exp [ - # - ^\ ^Hf^j} = 

^ k=l i,j ' 

iVexp { - £ r£>**** + (f^rgU^S?) ), (3.59) 



where 



.. 1 , 



fc=i 



I r 



r (2) 
"ll 



.(1) 
11 

I 



N = det A r = det (I - rffrg 1 ) 



(3.60) 



where I is the s x s unit matrix. 

The quadratic function in (|3.59)) can be represent in a compact form as: 



-V 2 A 12) &m 



with the (gi + q-i) x (qi + q 2 ) matrix r^ 12 ^: 

/ r (1) -v^lA- 1 ] r (1) -r^/l" 1 ! r (2] 

l 22 i 21 L^V J21 1 12 L 21 L^V J22 1 



,(12) 
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\ z 21 [ A r Jll 1 12 z 22 r 21 l^r J 12 1 
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and the inverse matrix A r 1 having a following form 



. 1 



\ [A- 1 )21 Ar ^22 



1 r ll r ll 



11 



(3.61) 



(3.62) 



t r (2) (1) 

1 r ll r ll 



The results in (|3.62|) reproduce for the considered situation the rules of 
the generalized star-product, derived in [7]. The matrix rj- describes the 
entire scattering process over the 2— vertex network. 

Now, let us note that the replacement of the local r( fc )-matrices by their 
unitary counterparts is equivalent to the unitary transformation of the mega 
S-matrix. Here we are giving the proof of that statement. We introduce 
formal notation § £ 5(rW)*£ = E p y ; ;,;'ew e:c i^i^u' V$ > denoting by V £ , ^ £ 
(^i, ^>z) the fermionic variables connected with the external (internal) links 
and by S(r^) the global scattering matrix. We shall denote the action in 

(|332])by A{r^; ^ e ,^e, *x, *z)- 

Then, performing a simple transformation of the Grassman variables 

jV( r (fc)) e -**S(rW)* e = j MjM^'iW 1 ^ = (3.64) 
.g-^(» -(fe)T ;xi!Xi 1 XfX£) — 
where r^j- = r^j, x 



,(i)J2) 



t r (2)_(l)^- 

1 ~~ r ll r ll 



(3.63) 



D X iD X ie- A{r{k)T '^' Xx ^ £X£) = N(r^ T )e- XeS(r(k)T)xe : 



x = (3-65) 
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we are coming to the relations: N(r^ T ) = N(A k ^) and S(A k ^ T ) = S(A k ^) T . 
The following relation takes also place: i>(r( fc ) ) = S(r^)* where * denotes 
complex conjugation. 

If r( fc > f = rW~\ then S^' 1 ) = S(r^* T ) = S(r^)* T . As S(r«) is 
unitary, if A k ^ is unitary matrix for each vertex k (see |16j). then 

5(rW _1 ) = S(rW) -1 . (3.66) 

and NirW 1 ) = N(r^)*. 

4 Bubble Diagram 

Let us apply the technique developed above to the calculation of the bubble 
network-diagram presented in Fig. 7. According to the prescription presented 
in the previous section the corresponding two particle i2-matrix is: 

where N is defined from i?^oi = 1> as it follows from the expression of two 
particle R-matrix (|2.5|) . This gives 

iV = E(-l) P ^ (3) of; 7 ^ (3) ^ 7l . (4-68) 




Y 

Figure 8: Bubble diagram constructed from two three channel vertices. 

The calculations can be easily done with the Mathematica program. In 
order not to overload the text with heavy expressions, we will only present 
hereafter the result for scale invariant R-matrices as defined in [6] , [llj , [12] . 
We begin by stating the result which is particularly simple, namely: 

Two scale invariant -matrices with scaling parameters Q\ and Q2 
respectively, are defining, via the expression (|4.67p . a scale invariant 
matrix with scaling parameter Q\Q2- 
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More explicitly, according to papers [BJ [TTJ [12] and relation (|3.33p , a scale 
invariant S-matrix defines the following parameterisation for the R-matrix: 



Hi 
a 



— , ru 
a 

(l + g)g 
a 

l + Q 



1 + Q 



T22 



ri3 



(l + g)g 

3 

a 

i + e 



g 

a a 



(l + g)g Q 

r32 = , ?"33 = , 

a a a 

1 + q + q 2 



(4.69) 



Substituting expressions (j4.69j) into (|4.67p and after renormalization of 
the product of the two i?( 3 )-matrices by N one will obtain: 



(2) 
11 



r (2) 
r 21 



2glg2 

1 + q\q\ 

-i + q\q\ 



,(2) 



1 



q\q\ 



12 



r (2) 
r 22 



1 + q\q\ 

2glg2 

1 + q\q\ 



(4.70) 



which is precisely the parameterisation of the scale invariant two particle 
S-matrix with scaling parameter q\ Q2 defined in (6j fTTT [T2] ! 

The bubble diagram we are considering is obviously the simplest non- 
trivial 2-vertex loop diagram. The 2x2 matrix r« can be constructed 

directly by use of both (|3.62j) and (14.670 formulas. Let us denote by 
(k = 1, 2) the two 3x3 matrices, each corresponding to one vertex. The 
numeration is taken so that the second channel of the first vertex (internal 
links) corresponds to the first channel of the second vertex, and the third 
channel of the first vertex is connected with the third channel of the second 
vertex. The first channel of the first vertex and the second channel of the 
second vertex correspond to the external links. The Ty matrix of the final 
two channel process can be expressed by the following expression [7]: 
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.(2) m 



21 



23 



A' 1 
A \2 



/J2) 
' r 12 
J2) 
'32 



V o 



21 



,(4.71) 



where 



.4 



12 



/ 1 



r (2) 
r ll 

Vr (2) 
V 31 



,(2) 
r 13 
r (2) 
r 33 



r (l) 
'22 

r (l) 
'32 

1 







1 J 



(4.72) 
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The normalization factor N equals to 



N 



+ 



Det[A 



12 



Det 




„(1U2) (1)(2) Jl) (2) 



22 ' 11 



„(1U2) 



(i)„(i) JiUm /U 2 U 2 ) 



' 22 r 33 



)( 



r ll r 33 



„(2)J2) 



„(2)^ 
13 

33 / 



(4.73) 



Note that the same expression (|4.73p for N can be obtained directly from 
formula (I4.68|) by use of Mathematica program. 



5 Triangle Vertex 

In the most general case, the three particle R-matrix is anisotropic and 
its elements depend on the channels 1,2,3. In such a case, it is necessary 
to specify, in the network diagrams, the channels of all the three particle 
vertices and follow which channels of neighboring vertices are connected. It 
is clear that all 3 3 = 27 types of connections in the triangle are in general 
possible and this will be defined by the considered problem itself. Here, 
as an example, we will present the case corresponding of the disposition of 
vertices as depicted in Fig. 8. Namely, the second channel of the first vertex 
is connected with the first channel of second vertex, the third channel of 
the second vertex is connected with the second channel of the third vertex 
and the first channel of third vertex is connected with the third channel of 
the first vertex. Notice that, according to Fig.5, the pairs of indexes (a,/3), 
(/?', a') and (7, 7') at the first, second and third position of R are correspond 
to the first, second and third channel respectively. 

Then, according to the prescription defined above, the i?-matrix corre- 
sponding to triangle diagram reads as: 

d(A)^Ti7 _ n(Ti5i5 3 ra/^fo 0717^3 ( _-\\P* A (r 74 \ 

The normalization factor iV is defined by the condition 

R (A) Z = I- (5-75) 

The matrix r^ 2 ^ , i, j = 1, 2, 3 describing the scattering over the triangle- 
graph in accordance with the formulas (|3.37p and (I3.55|) is: 

{ r$ \ 

rg> 



r (123) 



.(2) 



+ (5.76) 
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r 13 
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where r\- , k = 1,2,3 are the R-matrix parameters in the three vertices of 
the triangle diagram on Fig. 9. and 
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r ll 
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(5.77) 



with 



N 



det A 123 . 



(5.78) 



Again, the calculations of p^aa\ can be done with a simple Mathematica 
program. We will not present here the 20 big formulas, but just show the 
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normalization factor N: 



N 



+ 



E 
fi 



E>ci<5i<53 pi pO 1 CT3 / 1 \p(<T3) 

-"■0 1 o- 3 -"l a 1 S 1 n a 2 S 2 5 3 \ L > 
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'22 '11 
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(5.79) 



)JX) JX)JX)\ (J?) J® J2W1)\ / (3) (3) (3) (3) 



22 ' 33 



23 ' 32 



11 '33 



13 ' 31 



) H? 



22 



12 ' 21 



6 Conclusions 

We have presented a new approach to the calculation of complicated mega- 
S-matrices of the arbitrary network based on two ways: 

• Introducing R-matrices which are defined by the S-matrices of the local 
scatterings and taking the product of such R-matrices; 

• Formulating an action for the corresponding network by use of fermionic 
fields. 

In the examples of bubble and triangle diagrams we show the equivalence of 
this technique of computations with the generalized star-product procedure. 
The general proof valid in the case of an arbitrary graph was given in [16] . 

The action formulation has several advantages. First of all, it allows one 
to investigate network models with large number (N — > oo) and random 
disposition of vertices and formulate corresponding quantum field theory, 
which will boost further investigations. In this language one can take also 
sum over random dispositions of scatterers and analyze the problem of An- 
derson localization. Secondly, the approach makes the calculation of mega- 
S'-matrix elements of the regular, i.e. having some translational symmetry, 
networks very simple. Then the problem simply reduces to Fourier analysis 
of the action on a particular lattice defined by the network. This statement 
was demonstrated in Section 2 on a model of a chain of equidistant impu- 
rities. The third, and presumably most important advantage of the action 
formulation, is a possibility to generalize networks and consider models with 
multiparticle interactions. 
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